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Abstract: Gravitational scattering in the ADD-model is considered at both sub- and 
transplanckian energies using a common formalism. By keeping a physical cut-off in the 
KK tower associated with virtual KK exchange, such as the cut-off implied from a finite 
brane width, troublesome divergences are removed from the calculations in both energy 
ranges. The scattering behavior depends on three different energy scales: the fundamental 
Planck mass, the collision energy and the inverse brane width. The result for energies low 
compared to the effective cut-off (inverse brane width) is a contact-like interaction. At 
high energies the gravitational scattering associated with the extra dimensional version of 
Newton's law is recovered. 
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1. Introduction 

The ADD-model [1-3] is an attempt to solve the hierarchy problem, by introducing extra 
dimensions in which only gravity is allowed to propagate. For distances smaller than the 
assumed compactification radius, R , the gravitational potential will then be altered and 
has the form 



V(r) _ n r(n) . . 

-G NW H bn—j (1.1) 



m\m2 



where n is the number of extra dimensions, G^u\ denotes the ordinary 3+1-dimensional 
Newton's constant, S n = 27r ra//2 /r(n/2) is the surface of a unit sphere in n dimensions and 



1 Here we use the notations of [4], such that R is the radius and not the circumference. 
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r(n) is the Euler Gamma function. This implies that the strength of gravity increases 
much faster with smaller distance as compared with the normal 1/r behavior, and the 
fundamental Planck scale (related to the mass scale where the corresponding de Broglie 
wave length equals the black hole radius) is reduced and given by 

M D = (1.2) 

The presence of strong gravity at distances smaller than the compactification radius 
opens up for the possibility of observing gravitational scattering and black hole production 
at collider experiments and in cosmic rays. To eliminate the hierarchy problem, and not 
only reduce it, the new Planck scale should be of the order TeV, and LHC will be a quantum 
gravity probing machine. 

In order to quantify the amount of gravitational interaction, the theory was formulated 
as a field theory in [5,6]. As the extra dimensions are compactified, the allowed wave num- 
bers (and hence momenta) in these dimensions are quantized Kaluza-Klein (KK) modes. 
The KK modes can of course enter both as external and internal particles in the Feynman 
diagrams derived from the theory. When the KK modes are internal (as for elastic gravi- 
tational scattering) they have to be summed over. The problem is that the sum over KK 
modes diverges for 2 or more extra dimensions, 

Y.z^T¥ = s ^h^TB dm - (L3) 

I 1 

Here I enumerates the allowed KK modes with momenta mj in the extra dimensions, 
m = \mj\, and k is the exchanged 4- momentum in our normal space. (We will for simplicity 
call this object a propagator, despite the fact that the Lorentz structure is not included.) 

In the original papers [5, 6] this divergence problem was dealt with by introducing a 
sharp cut-off, M s , argued to be of the same order of magnitude as the Planck mass, as new 
physics anyhow is expected to occur at the Planck scale. Various mathematical forms of 
cut-offs have also been discussed in [7]. For n > 3 and momentum transfers small compared 
to M s , the sum was then estimated to give 

1 1 M n ~ 2 

R n M™~ 2 « (1.4) 

n-2 s G N{4) {n-2) Ml +2 V ; 

In the Born approximation this would lead to a cross section of the form [8] 

da s 3 / M n ~ 2 \ 2 

- — rrr ~~~ryj Hspin, z) (1.5) 

dz {n-2) 2 \Ml +2 ) K 1 ' V ' 

where z is cosine of the scattering angle in the center of mass system, ^fs the total cms 
energy, and F a function taking spin dependence into account. 

Ordinary gravitational scattering in 3 + n dimensions would correspond to a poten- 
tial oc l/r( n+1 ), but the scattering given by eq. (|1.5| ) has a completely different angular 
behavior. In particular the expected forward peak is totally absent. Fourier transforming 
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the amplitude in eq. (|1.4j ) to position gives a (S-function potential, oc 5(f), and the corre- 
sponding Born approximation cross section in eq. ( p-5| ) is therefore isotropic. Thus it is 
obvious that the approximation in eq. ( |1.4| ) does not contain the full story of gravitational 
scattering in the ADD-model. 

An attempt to solve this problem has been presented by Giudice, Rattazzi, and Wells 
[4]. These authors point out two important facts: 

i) For an interaction with a large Born amplitude but a short range, the cross section 
is not determined by the Born term alone. Higher order loop corrections reduce the cross 
section and guarantee that the unitarity constraint is obeyed. 

ii) The constant term in eq. ( |1.4| ), which represents a dominant part of the amplitude 
in eq. ([L.3j), corresponds to a contribution to the cross section from zero impact parameter, 
and should therefore give a negligible contribution to the cross section, at least when the 
incoming wave packages do not overlap. Consequently the important part of the amplitude 
in eq. ( |1.3[ ) must in this case be the smaller fc-dependent terms, which have been neglected 
in eq. ([111) . 

In case the interaction is dominated by small angle scattering the cross section can be 
calculated in the eikonal approximation, in which the all-loop summation exponentiates 
[9—11]. The cross section is then given by 



a el = J d 2 b ± \(l-e 



fftot = J d 2 b ± 2Re(l-e i ^) (1.6) 
with X (b ± ) = 1| ^c^AbUqI)- (1-7) 

Thus, if the absolute value of the eikonal, x, is small compared to 1, we in general expect 
small corrections from the higher order loop contributions, while for large x -va l ues the 
cross section saturates, and the effective integrand in eq. (|l.6|) is close to 1. We also note 
that when \ 1S real, the scattering is purely elastic. In this paper we will focus on elastic 
collisions mediated via (multiple) exchange in the t-channel. 

It is also pointed out in [4] that in the eikonal limit the Born amplitude does not 
depend on the spin of the colliding particles, and is therefore universal. Expressed in the 
fundamental Planck mass Mb in eq. ( |1.2| ) it is given by [4] 

In [4] a divergent part is subtracted from the integral in eq. ( |1.3| ) or ( |1.8| ) using di- 
mensional regularization. This subtracted part corresponds to a narrow potential localized 
at f = 0. Although the remainder is singular for n equal to an even integer, its Fourier 



transform (the eikonal % m eo i- (hi)) 1S finhe everywhere. Assuming the eikonal approx- 
imation to be applicable in the transplanckian region s 3> M D , the authors of [4] thus 
obtains a reasonable result, where the gravitational scattering cross section grows with 
energy cx (s/M" +2 ) 2 / n . However, we ought to be worried by the fact that the part of the 
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amplitude, remaining after the subtraction, grows for larger momentum transfers, and is 
largest for backward scattering. This implies that the conditions for the eikonal approxi- 
mation are not satisfied. The formal problems with divergent integrals also indicate that 
this result could be regarded as based more on physical intuition than on a solid theoretical 
foundation. These uncertainties also make it difficult to estimate the limit beyond which 
the result should be applicable, and how the gravitational scattering behaves for lower 
energies. 

In this paper we want to study in more detail the result of various physical effects, which 
can tame the divergences. These effects give effective cut-offs for high-mass KK modes at 
some scale (here referred to as M s ), which does not have to be the same as the Planck 
scale Md- Our result does indeed confirm the relevance of the eikonal approximation and 
the result in [4] at very high energies. For lower energies the behavior is different, wide 
angle scattering is dominant and the amplitude does not exponentiate. Instead the all-loop 
summation gives a geometric series. This implies that there will be a change in the energy 
dependence, and for lower energies the cross section varies more rapidly, proportional to 
ocs 2 M* n - 2 /M 2 D n+ \ 

We want to emphasize that in this paper we do not discuss phenomena like black hole 
formation or other nonlinear gravitational effects, which are expected to modify the final 
states for very high energies and central collisions. For a discussion of such effects we 
refer to ref. [3,4,12-17]. We also neglect possible interference with strong and electro- weak 
effects and we study reactions for non-identical particles such that KK modes appears only 
in the t-channel. Some remarks on s- and u-channels are however made in sees. 3.3-|6.5. 

The approach in [4] will be discussed in more detail in sec. In sec. ||, we will introduce 
a finite width of the brane, on which the standard model particles are assumed to live, and 
see how this leads to a finite amplitude. A similar effect is obtained by assuming that 
the position of the brane is not fixed in the extra dimensions [18, 19]. Fluctuations in 
the brane then result in a kind of surface tension or "brane tension". The Born term 
is discussed in sec. |2| and higher order loop corrections in sec. ||. Here we also study in 
which kinematical regions the Born term dominates, where the eikonal approximation is 
valid, and the behavior of the cross section in regions where the scattering is approximately 
isotropic. The results for scattering cross sections in those different kinematical regions are 
then presented in sec. §. Finally we will summarize and conclude in sec. 0. 



2. Problems and divergences 



The integral in eq. (1.3) or ( |l.8|) is divergent for n > 2 and n < 0, but converges for 
n- values in the intermediate range < n < 2. To give a physical meaning to the integral 
for n > 2, a finite result can be obtained by analytic continuation from smaller n-values, 
corresponding to a dimensional regularization. The resulting amplitude, presented in [4], 
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is given by the expression 2 



-Trtrril-- 
2 



M D 



Ml 



(2.1) 



We see that this expression is finite for odd integers n, but singular for even re, where the 
T-function has poles. 

The result in eq. (|2.1| ) is equivalent to a subtraction of terms, which are proportional 
to ^-functions or derivatives of (^-functions at f = 0, and therefore may be expected to give 
negligible contributions to the cross section. Inserting eq. (^^) into the two-dimensional 



Fourier transform in eq. (1.7), we see that this integral is also divergent. It can be given a 
finite result by introduction of a convergence factor: 



X 



with b r 



s(47r)f- 1 r(n/2) 



2M£+ 2 



l/n 



(2.2) 



We note that although the amplitude A^om in eq. (^j]) is singular for even re, \ is finite. 
Thus x(b) (like the potential V(r), to be discussed below) can be analytically continued 
to finite values for all n-values. (A finite amplitude, which corresponds to a potential 
proportional to \/r n+l for n even, is oc (— k 2 ) n l 2 ~ l In (— k 2 ).) 

The result in eq. ( |2.2| ) is a single power oc 1 jb n , and the scale factor (or characteristic 
impact parameter) b c is defined so that M = 1 when b = b c . If this expression is inserted 



into eq. (1.6), we see that the term quadratic in x, which is the Born term, dominates 
the integrand for b > b c , where x < 1> but higher order corrections are important in 
constraining the scattering probability for b < b c . 



In ref. [4] it is assumed that eqs. (1.6) and (2.2) should give a realistic approxima- 
tion to gravitational scattering in the transplanckian region s ^> (apart from special 
effects like black hole formation, which are treated separately). The net result is then 
that the total scattering cross section grows with energy proportional to b 2 , or equivalently 
oc (s/M" +2 ) 2 / n , (cf. eq. Q below). 

The exponentiation in the eikonal approximation in eq. ( |1.6[ ) follows when the scatter- 
ing is dominated by small angles [9-11]. The one-loop contribution is then dominated by 
its imaginary part, and the all- loop summation gives an exponential. 

The one-loop diagram in fig. |l|a is given by the following expression: 



A\— loop(^ 



(2tt)< 



■-4Born(g 2 )^Born((A; 



(pi - q) 2 (p 2 + q) 



2 ■ 



(2.3) 



Here p± and p2 denote the momenta of the incoming particles, the total momentum ex- 
change is k and the loop momentum q. 

The imaginary part of the integral in eq. ( |2.3| ) is coming from on-shell intermediate 
states (denoted i in fig. |l]), and can be calculated using the Cutcosky cutting rules. This 



2 We have here inserted a minus sign not present in [4] 
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Figure 1: (a) The one-loop contribution corresponding to exchange of two KK modes. The 
KK modes are drawn as thick lines and standard model particles as thin lines, (b) The two-loop 
contribution. 



implies that the two propagators in eq. (^^) are replaced by 5- functions, which (with the 
approximation q 2 « q\) gives the result 

/d 2 ~ 
j^A Born (-ql)A Born (-(k ± - q ± ) 2 ). (2.4) 

If ^Bom falls off for large k±, the Fourier transform to impact parameter space of the 
one-loop contribution is proportional to \ 2 {b±). The sum over multi-loop ladder diagrams 
with different number of KK exchanges exponentiates to (ix — X 2 /2 + •••) = e* x — 1, and 
the all order eikonal amplitude is given by 

A cik (k 2 ) = -2is J d 2 b ± e fk ^(e ix - 1). (2.5) 

With the Born amplitude in eq. fl2.1| ) we have, however, some problems. First the real 
part of the integral in eq. (^^) is not small and negligible compared to the imaginary part. 
It is strongly divergent for n > 2, as A Born increases proportional to q n ~ 2 or q n ~ 2 ln(—q 2 ) 
for large q. Secondly the integral in eq. ( |2.4[) should only go over physical intermediate 
states, which means q± < yfs/2. The Fourier transform of the convolution in eq. ( |2.4| ) 
gives x 2 only if the integrand in eq. (|2.4j ) falls off so rapidly, that the integral effectively 
can be extended to all q±. This is not the case with the amplitude in eq. (|2.1| ). 

We conclude that, although the result in eq. (|2.2| ) and (|1.6| ) is an intuitively reasonable 
result for scattering in a rapidly falling potential, it should be worrying that it is derived 
from an amplitude, which grows for large momentum transfers and large scattering angles, 
while the eikonal approximation is proven to be valid only when scattering at small angles 
is dominating. 

At the root of this problem lies the fact that the subtraction, which gives the amplitude 
in eq. (|2.l| ) and is a result of the analytic continuation in the number of extra dimensions, 
does not automatically remove all parts corresponding to 5- functions at f = 0. The def- 
inition of the potential as the Fourier transform of eq. ( [2.1| ) is problematic. To illustrate 
this we study the most simple example represented by the case n = 3. In the rest frame 
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we have ko = and k 2 = —k 2 . The integral in eq. ( |1.8| ) is then proportional to 

m? dm f {m 2 + Ar — k 2 )dm f - 2 f dm 

k 2 + m 2 J k 2 + m 2 J J k 2 + m 2 

The first term, the integral, represents an infinite subtraction. Its three-dimensional Fourier 
transform gives a 5-function at r = with an infinite weight. The second term corresponds 



to the result in eq. (2.1). We may try to define its Fourier transform V{r) as a distribu- 
tion in the standard way, multiplying with a test function and interchanging the order of 
integration. For test functions of the form exp(— ar 2 ) we then get (with k = \k\ and the 
constant C appearing in eq. Q4.5|)) 



7T , 



d 6 re~ ar V{r) = C{ ) / d 6 kk / d 6 re 



-ar 2 irk 



= -C* d*kke- k2 ^£)l 
2 J V 

= -CX^^fd. (2.7) 

We note that this result is finite, and goes towards when the test function approaches 
a constant, i.e. when a — ► 0. For r / we find V(r) = C4-7T 2 /r 4 by Fourier transforming 
from k to f using a convergence factor. Integrating this contribution with the test function 
above, we get the divergent result 

C167T 3 r e ~ ar2 %. (2.8) 



Thus this definition, V(r) = C4tt jr for r > 0, is incomplete since the result in eq. Q2.7| ) is 
finite while the integral in eq. ( |2.8| ) is infinite. It looks as if a ^-function, 6(f), with infinite 
weight is missing. 

We conclude that the separation in eq. (|2.6| ) does not in itself remove all terms related 
to 5- functions at f = 0. Instead we argue in the next section that dynamical effects will 
remove the divergencies and give finite results. 

In the next section we will consider possible mechanisms, which can suppress high KK 
masses and give an effective cut-off to the integral in eq. (|1.8| ) . These mechanisms have real 
dynamical motivations, and we will see that such finite cut-offs do remove all divergences 
and give well-defined results. For high energies the Born amplitude indeed falls off for large 
momentum transfers, and the eikonal approximation is applicable. For lower energies this 
is not the case, and we will in sees. || - ^ discuss the resulting amplitudes and cross sections 
for different relations between the energy, the Planck mass, and the cut-off scale. 

3. Possible solutions 

In the ADD-model the standard model particles are assumed to live on a thin brane. The 
mechanism behind this assumption could possibly be taken from string theory [3], but is 
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not a part of the ADD-model itself. The problems discussed in the previous section are 
related to contributions from KK modes with very high masses. In a relativistic quantized 
theory there are also formal problems with an infinitely thin and infinitely rigid brane. If 
the brane is not infinitely thin, but has a finite width, this will effectively suppress the 
coupling to high-mass KK modes, with wavelengths shorter than the brane width. If the 
brane really is infinitely thin, then it must be impossible to determine its position with 
infinite accuracy. In [18, 19] it is demonstrated that the fluctuations in the position of 
the brane suppresses high-mass KK modes, in a way similar to the effect of a finite brane 
width. The emission or absorption of a KK mode gives a recoil to the brane, and the 
fluctuations in the location of the brane can then be regarded as a result of an effective 
"surface tension" in the brane. 

Let us for definiteness study the effect of the assumption that the standard model 
fields penetrate a finite distance into the extra dimensions, which gives an effective finite 
width to the brane [20]. (The possibility of fluctuating branes, studied in [18,19], give 
similar results, albeit with a different physical interpretation.) To be specific we assume 
a Gaussian extension, but this assumption is not essential for our conclusions. Thus we 
assume the standard model fields to have a wave function with the extension 

— 

M s \ 2 __ja M a /4 



m=[^= e-r^t* (3.1) 



into the extra dimensions, with y denoting the coordinate in the extra dimensions. The 
overlap between two standard model fields and a KK mode of mass m (what we have in a 
vertex) is then proportional to 



/ 



Z 2 



d - e m-y ( MA' * e -y*Ml/2 = e -«W) (3.2) 
2tt ) 



or, in other words, the squared absolute value of the wave function in y-space Fourier 
transformed to m-space. The exchange of a KK mode will have this suppression factor 
occurring twice, once at every vertex. In total the exchange of a KK mode with mass m 
will therefore contribute to the sum in eq. (fbH) with a suppression factor 



e -™ 2 /M?_ (3.3) 

Implementing the physical requirement that the standard model particles live on a 
narrow brane does therefore in itself imply a finite "effective" propagator, 

irs n [ |p^ e -mv^ (3 . 4) 

J k z — m z 

for the exchange of 4-momentum k. (The factor R n comes from the density of KK modes 
and S n = 2-7r n//2 /r(n/2) is the unit surface of a sphere in n dimensions.) We note in 
particular that this expression (in contrast to the expression in eq. ( |2.1[) ) falls off like 1/k 2 
for large momentum transfers, such that —k 2 ^$> M 2 . This implies that for high energies, 
s S> M 2 , t-channel interaction is dominated by small values of —k 2 /s, i.e. by small angles. 
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In the following sections we will show that the Fourier transform of the propagator 
in eq. ( |3,4|) gives a potential, which falls off oc l/r n+1 for distances larger than the brane 
width, given by 1/M S , and smaller than the compactification radius. Outside this range, 
both for r < X/M s and for r > 2nR (where the massless graviton dominates), it varies 
oc 1/r. We will also study the resulting scattering cross sections under different kinematic 
conditions. 

4. The Born term 
4.1 Amplitude 

As described in section ^, several physical mechanisms result in effective cut-offs for high 
masses in the KK propagator. After multiplying eq. (3.4) by the coupling 4-7rG7v(4), con- 



tracting Lorentz indices (not explicitly included here), and using the relation G N ,^ = 
we get the following result for the Born amplitude for ultra-relativistic small 
angle scattering: 



W*) = T^S„r^f^e-V«l. (4.!, 



M£+ 2 Jo k 2 - m 2 

For large angels there are less important corrections from spin polarization which we neglect 
here and in the following. This integral is convergent and finite for all negative values of 
k 2 = t (including when n > 3). It is easy to find the result in the limits of large and 
small (negative) t-values. 

• Large momentum transfers; —t ^> M 2 

When —t is large compared to M 2 , the term m 2 in the denominator in eq. (4.1) can 



be neglected, which gives the result: 

Thus for large momentum transfers (larger than the cut-off) the Born amplitude falls 
off proportional to 1/t. 

Small momentum transfers; —t <C M 2 

For smaller t, and n > 2, the integral is dominated by m- values of the order of M s , and 
therefore t can now be neglected in the denominator. We then get the approximately 
constant result: 



A 



Born 



(*) 



M 



n+2 



dmm n - 3 e- m2 / M ° 



2W 2 / AL 



(n-2) \M D 



M 2 D M 2 ' 



(4.3) 
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Thus for momentum transfers, which are small compared to the cut-off, the Born 
amplitude is approximately constant for n > 2. For n = 2 the result for small t has 
instead a slowly varying logarithmic dependence, proportional to ln(— M 2 /t). 

4.2 Potential 

To get the classical non-relativistic potential we start directly from the effective propagator 
in eq. fl3.4|) multiplied with the coupling constant 47rCrjv(4) • Going to the rest frame, where 
ko = and t = —k 2 we find the corresponding potential as the three-dimensional Fourier 
transform: 



V(r) 1 f d 3 k 



/d k jfoy / T2\ 

(2^)3 e ABmn{ -~ k > 



m\m,2 2s 2 

i. Or, / , „_i _ rr ,2n\,f2 / U hot. 



dmm n - l e- m ' M * 



2M™ +2 (2^) 3 7o J m 2 + k 2 



1 Q /*°° £>~ 



-mr 



^ ^2vr 2 / dmm"-^'.- . (4.4) 

2M™ +2 (2^) 3 y r V ; 

This represents a weighted sum of Yukawa potentials. The integral can be expressed in 
terms of error functions, but we are here primarily interested in the behavior for large and 
small values of r. 

• Large distances; r > 1/M S 

For distances larger than the brane thickness the integral is effectively cut off by the 
factor e~ mr , and the result becomes insensitive to the Gaussian cut-off e _m2//A/ ». It 
is then approximated by 

V(r) _ -1 S n f°° n _! e~ mr S n T(n) 1 

dmm . __ = . (4.5) 



m x m 2 2M™ +2 4vr J r 8vrM™ +2 

We see that for distances large compared to the brane thickness (but small compared 
to the compactification radius) we recover the result from eq. Ql.l| ), a potential falling 
off proportional to l/r n , corresponding to the expected (3 + n)-dimensional version 
of Newton's law. When r is increased, smaller m- values ~ 1/r are important in the 



integral in eq. (4.4) or (|4.5|). The phase space factor m n then gives this power-like 



fall off for distances large compared to M s . 
Short distances; r < 1/M S 

For smaller r-values we find instead that the factor e~ mr is irrelevant, and the result 
is 

nr) _ yr^,,,--,.!^*.! (4 . 6) 



mim 2 2Mp +2 47r r J 87r M^ +2 r 
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Due to the cut-off, the integral in eq. (fi.4| ) is dominated by m-values close to 
M s for all r-values smaller than X/M s . Thus, when the distance is smaller than 
the brane width, the result is a potential proportional to 1/r, corresponding to 
a standard 3-dimensional Coulomb potential, although with a coupling constant 
~ M"/M™ +2 ~ (M s R) n G N{4:) instead of G N{A) . Thus the coupling is enhanced 
by a factor ~ (M a R)» = ( ^P^ification ^radius )n 

4.3 Eikonal 

In a similar way we can calculate the eikonal x(b) by a two-dimensional Fourier transform 
in the transverse coordinates: 



S Sn r dmm n - 1 e- m2 / M » [ d 2 k ± e fk ^- 1 



2M n D +2 (2vr)2 J J w ^ m i + kl~ 

-s S n f°° x _ m 2 /M 2 f°° kj_dk± 

7 2tt / drum e ' s / — 5 -s- Jo(k±b) 



2M£ +2 (2vr) 2 y y m 2 + A; 2 

C /"OO 

dmm n - l e- m2 / M sK Q {mb). (4.7) 



2M£+ 2 2tt 7 

This integral can be expressed in terms of confluent hypergeometric functions of the 
second kind: 



X(b) 



sM? n 7T 
M™+ 2 [ 2> 



n/2-l 



2a2 



l 2' ' 4 



(4.8) 



This expression can easily be used in numerical calculations. For an intuitive picture, the 
result for large and small 6-values can be estimated in the same way as the approximations 
in eqs. (4.5, 4.6). 



Large impact parameters; b S> 1/M S 

For large arguments the asymptotic behavior of the Bessel function Ko(mb) is pro- 
portional to exp(— mb) / V mb. This implies that for large b the Gaussian cut-off is 
unessential, and we find the eikonal 



X Q>) « * [ dmm n - l K (mb) = ^^2™~ 4 T 2 (n/2) • 1 

M ' 2M™+ 2 2tt J UV ' M™+ 2 vr y 1 ' b n 

(4*)*"* - •-. (4. 



Small impact parameters; b -C 1/M S 

For small arguments we have Ko(mb) ~ ln(l/(m6)), which implies 
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X(b) 



2M™+ 2 2vr 



/ 



7T2 



4 M 2 D \M D J V 



where VKf) is the psi or digamma function. 



1 

mb 



(ln(M s b) + 



(4.10) 



Thus we see that the eikonal falls off oc l/b n for large b, and grows logarithmically 
when 6^0. Using the quantity b c from eq. (^^) and keeping only the dominant term 
ln(M,6) in eq. (4.10), we can write the results in the form 



X(b) 

X(b) 
with bd 

and b r 



i \ n 



2 n " 1 r(n/2) \M s b 
1 
Ml 



; b<b d 



s(47r)5-T(n/2) 



2M 



n+2 



1/n 



(4.11) 
(4.12) 
(4.13) 

(4.14) 



The separation line 6^ = 1/M S is an estimate of the b- value where x(b) changes behavior. 
As an example fig. shows these approximations for x together with the exact result for 



3 and -^/i = Md = 1 in units such that M s = 1. As x is proportional to s/M^y , a 
change in s and/or just corresponds to a shift of all curves the same distance up or 
down. 



5. Higher order loop corrections 

We note that three different energy scales enter the expressions for the Born amplitude 
in eqs. ( fO| , |4.8|) : ^/s, Md, and M s . Here yfs is the total energy in the scattering, Md 
is the fundamental Planck scale determined by the compactification radius R, and M s is 
related to the width of the brane (or the brane tension). The result depends on the relative 
magnitude of these quantities, and in the following we will successively discuss five different 
kinematical regions, which are illustrated in fig. [3|. 

5.1 Eikonal regions, s ^> M 2 

We study the scattering process in the overall cm system, where the momentum exchange 
has no 0-component, k = (0;k) and t = —k 2 . From eq. Q4.2] ) we see that ^Born falls off 
oc l/£r for k 2 > M 2 . Thus for high energies, such that s M 2 , corresponding to region 1 
and 2 in fig. ||, the Born term is dominated by small values of k 2 /s, i.e. small angles. This 
implies that the eikonal approximation is applicable. We note in particular that it is M s 
rather than Md, that sets the scale for when the eikonal approximation is relevant. The 
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ln(-x) 




Figure 2: The logarithm of \x\ as a function of impact parameter for n = 3. The curves show 
the example where y/s, Md, and M s have the same magnitude, and the units are chosen such that 
y/~s = AId = M s = 1- This also implies that 6<j = 1. The uppermost line is the large b limit of x 



taken from eq. ( 4.11 ) and the lowermost line is the small b limit of x taken from eq. ( 1.12 ). The 



interpolating line is the exact expression eq. (4.8). Note that a change in s and/or Mu, keeping 



M s constant, just corresponds to shifting all curves up or down. 

one-loop contribution is here dominated by its imaginary part, obtained when the particles 
in the intermediate state i in fig. [l]a are on shell. The contributions from multi-loop ladder 
diagrams (fig. [j]b) exponentiate, and the scattering amplitude is given by eq. Q2.5] ): 

A cik (k 2 ) = 2is J d 2 b ± e fk ^(l - e ix ). (5.1) 

From eq. ( |5.lD we see that the higher order corrections are important when x ls OI 
order 1 or larger. Correspondingly the Born term dominates when \x\ < 1- We see from 
eqs. ( gl| pUD that X varies only logarithmically when b is decreased below the point 



b = bd- The importance of the higher order corrections for the integrated cross section 
therefore depends on whether or not > 1. This relation is satisfied whenever b c > b^, 

or equivalently when s > s c d, with s c d given by 

2 M n+2 

Scd= (An)^r(^) ' ^ 

This defines the boundary between region 1 and region 2 in fig. |3L In region 1 higher order 
terms are important for b < b c , and the exponentiation in eq. ([O]) is essential to keep the 
amplitude within the unitarity constraints. 

The difference between regions 1 and 2 is illustrated in fig. |j. Fig. |4|a corresponds to 
region 1, where the energy is high, and b c > bd- The absolute value of the eikonal x is 
smaller than 1 for b > b c , and in this range the approximation in eq. ( [4.1 1| ) is relevant. For 
b < b c , \x\ is large and rapidly varying, which causes the exponent in eq. (|5.1| ) to oscillate 
rapidly. 
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Figure 3: The (y/s, M s )-plane for n = 3 and Mb = 1. The straight line separating region 1 and 
4 is yfs = M s while straight line separating region 4 and 5 is the line where the real and imaginary 



parts in eq. (5.4) have equal magnitude. The power-like curve separating region 1 and 2 is s c d from 
eq. ( |5.2| ) as a function of M s and the line separating the regions 4 and 5 from region 3 is the line 
where I^Bom-X"! = b sec cq. (5.4). In the regions 1 and 2 s is larger than M s , and, at least for 
s 3> M s , the eikonal approximation is correct. In region 1 the eikonal is, depending on b, either 



large compared to 1 or given by eq. ( 4.11 ). In region 2 on the other hand the 6-range where x lS 
small includes a region where it is described by cq. ( 4.12j ) . In region 3 the correction corresponding 
to higher order loops is small, but in region 4 it is important to assure unitarity. 



Fig. 0b corresponds to region 2. Here |x| < 1 except in a very small region 



b < "tt exp 



4M"+V- 



(5.3) 



around the origin. Therefore the Born term dominates the cross section, and higher order 
terms give only small corrections. 

5.2 Non-eikonal regions, s < M 2 



The Born amplitude in eq. (|4.3| ) is almost independent of the momentum exchange k when 
k -C M s . When \fs -C M s (regions 3, 4, and 5 in fig. ||) this includes all kinematically 
allowed A:-values, which implies that the scattering is almost isotropic. Thus the exchange of 
the KK modes corresponds effectively to a contact interaction. (For wide-angle scattering 
we also expect corrections from spin polarization. This effect is neglected in the following.) 
The one-loop contribution in fig. [l]a is then represented by the diagram in fig. ||a, which 
is easily calculated. We denote the momenta in the intermediate state P/2 ± q, with 
P = p\ + P2, as indicated in fig. ||a, and in the cms we have P = (W, 0). The vertices are 
then given by the Born term in eq. ( |4.3D , with an effective cut-off when q = M s . The result 
is therefore 



loop(fe 



/ 



d A q 
(27T 



-A 2 

4 Born 



(P/2 



(P/2 + 
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Figure 4: — \ as a function of impact parameter for two examples with n = 3. (a) High energies 
corresponding to region 1 in fig. [| with y/s = 10 TeV, M s — 1 TcV and Mn = 1 TeV. The upper 
curve is the approximate expression in cq. ( |4.9| ), and the lower curve the exact expression eq. (5.4). 
(b) Kinematics corresponding to region 2 in fig. H y/s = 0.1 TcV M a = 1 TcV and M D = 1 TcV. 



The upper curve is the approximate high b expression in eq. (4.11), the lower curve the approximate 



low b expression in cq. (4.12) and the interpolating line is the exact expression in eq. (4.8) 




Figure 5: When the exchanged momentum is small compared to M s , the KK propagators are 
effectively replaced by vertex factors. The diagrams in fig. |l| can then be drawn as above with only 
standard model particle lines. 



1 AM 

= Al oin -X, withX«— ^(ln s - + iv) (5.4) 

We note here in particular that the result is a constant, independent of the momentum 
transfer k. Therefore also the one- loop amplitude can be effectively regarded as a contact 
term with a cut-off when k > M s . Consequently the two- loop diagram in fig. |B|b can be 
calculated in the same way as the one-loop diagram, and the result is 

^2-loop = ^l-loop • ^.Born^ = ^.Born ■ (^-Born^) 2 • (5-5) 

In the same way we can calculate ladder diagrams with more loops. Summing all ladders 
we obtain 

^-ladder = ^.Born (1 + ^Born^ + (^Born^) 2 +...) = T~~~~ Tp- (5.6) 

t — /iBorn^ 

We see that instead of the exponent in the eikonal regime (where forward scattering domi- 
nates) we here obtain a geometric series from ladder type contributions. The importance of 
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(b) 



Figure 6: (a) An example of a non-ladder diagram contributing to the elastic cross section in 
region 5 in fig. [s| (b) An example of a diagram contributing to the inelastic cross section in region 
5 in fig. |. 



higher order corrections is now determined by the quantity ^Bom^- When |ABom^1 ^ 1 
the Born term dominates. This corresponds to region 3 in fig. ||. 

When instead l-ABom-Xl > 1, we expect different results depending on whether it is 
the real or the imaginary part which dominates. When ln(4M^/s) < tt, the imaginary 



part dominates the loop integral in eq. (5.4). Thus this diagram is dominated by real 
intermediate states i in fig. |l]a, and the ladder diagrams in fig. [l|b or fig. [|b should be 
important higher order corrections. This corresponds to region 4 in fig. S 

When ln(4M^/s) > it (region 5 in fig. ^) the real part dominates the loop integral. 
This implies that inelastic scattering and virtual intermediate states are essential. We 
then expect important contributions from more complicated, non-ladder, diagrams, like 
the examples shown in fig. |6|. For this reason we do not expect the result in eq. ( |5.6| ) to 
be representative for a sum of all higher order corrections in this kinematical region. 



6. Cross sections 



Below we successively discuss the cross sections obtained in the five different regions in 
fig- | 

6.1 Region 1, s > M s 2 and |x(MI > 1 



In this region s > and xipd) > 1- As discussed in section 5.1 the scattering is suppressed 
for —t > Mg. The first constraint therefore means that the cross section is dominated by 
small angle scattering, the imaginary part dominates the onedoop contribution, and the 
eikonal exponentiates. The cross section is then given by 



a = / d 2 b ± 2Re(l - e 



(6.1) 



The effect of the constraint |x(6<i)| > 1 was illustrated in fig. ||a. It implies that < b c , and 
that the approximation x ~ —(b c /b) n in eq. ( 4.11 ) is relevant for all b > b c . In particular 



this means that, for b > b c (> bd), we have \x\ < 1 and 2Re(l — e* x ^ ± - ) ) ~ x 2 - For central 
collisions, where b < b c , higher order loop corrections are important to satisfy unitarity. 



- 16 - 



Here |x| is larger than 1 and rapidly varying, the exponent in eq. (6.1) is oscillating, and 
therefore 

' 2Re(l - e ix ^) \ « 2. (6.2) 



(For b < bd the variation in \ is only logarithmic, and therefore small in relative magnitude. 
As x is very large, the variation may still be large in absolute magnitude and such that 
eq. Q6.2p is valid also here.) Inserting these results into eq. (|6.1|), we get (for n > 2) the 
following result for the total cross section: 

cr« f bC d 2 b-2+ r d 2 b(^) 2n = TTb 2 c (2 + —) = 2Trb 2 c ^^. (6.3) 
Jo Jb c \bj n-Y n-1 

When s is increased, a grows proportional to b 2 oc s 2 l n . We note that the cross section 
is dominated by central collisions with b < b c (especially for large n), with only a small 
contribution from larger impact parameters. Integrating the constant 1 in the parenthesis 



in eq. (5.1) between and b c gives a dominant forward peak, with oscillations at larger 
angles. The amplitude for these oscillations falls off proportional to 1/k 3 ' 2 , corresponding 
to da/dt oc l/k 3 for the cross section. 



For large k the dominant contribution in eq. (5.1) comes from the term e tx and a small 
range of 6- values around b s , where 

f n \ "+ 1 



For these b- values the frequencies of the exponents e lk ' b and e %x< ^ in eq. (5.1) oscillate in 



phase, which gives an enhanced contribution. Using the saddle-point approximation we 
get from this contribution (apart from logarithmic corrections) a cross section which falls 

n + 2 

off like da/dt oc l/£ "+ 1 . This contribution is dominating for k > n/b c , where |x(6 s )| > 1- 
As pointed out in [4] it corresponds to classical scattering in a l/r n+1 -potential. For 
small scattering angles 8 we have for a non-relativistic particle with mass mi moving with 
constant speed v and momentum p = m\v in the potential of a mass mi 



= -L / dtF ± (r) (6.5) 



-Gjsr(4)R n S n r(n 



vri\m2 f°° dr d ( 1 x 1 1 



m iv J-oo v db V Vr 2 + b 2 
n(2V¥) re r(§) m 2 1 
8ttv 2 M% +2 b n+1 ' 



From this we see that if mi = s/(4m2) 



-*nonrel 



l 

' (6.6) 



4vpb, 

agreeing parametrically with eq. fl6.4|) . (A numerical difference is expected since eq. (fc 
is ultra-relativistic whereas eq. ( |6.6| ) is a non-relativistic result.) This behavior is discussed 
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in more detail in [4], and we note that in this region, where s is much larger than both 
M 2 and s c d, our result is consistent with the result of this reference. A necessary condition 
is, however, that y/s, Md, and M s have values such that b s > bd = 1/M S , which for fixed 
k- value gives a minimum value for M s . If this relation is not satisfied, the phase variation 
in exp(ix) is given by eq. ( |4.12| ) rather than eq. ( |4.11 ), and therefore we do not get the 
phase coherence in the integral in eq. (|5.1|). 



6.2 Region 2, s > M 2 and \x(b d ) 



< 1 



In region 2, s is larger than M 2 but smaller than s c d, and therefore bd > b c . A typical 
example is illustrated in fig. |4|b. We see here that \x\ is small compared to 1, apart from the 
logarithmic peak for very small b. The influence of the small b peak is also suppressed by 
a phase space factor proportional to b db. The cross section is therefore well approximated 
by the Born amplitude. 

The largest contributions to the cross section come from b- values in the neighborhood 
of bd] for larger b, x fells off oc (b c /b) n , and for smaller b the scattering is limited by the 
smaller phase space ~ bdb. These b- values are just in the transition region between the two 
asymptotic forms in eqs. ( [4.11 , 4.12). To get a good estimate of the cross section we should 



therefore use the exact expression for x i n e 1- Q4.8| ). For an order of magnitude estimate 
we may, however, approximate x by the asymptotic result x ~ —(b c /b) n for b > bd, and by 
a constant = — (b c /bd) n for all b < b d - This gives the following qualitative estimate for the 
total cross section: 



hd .0, (b r \ 2n r ,o. (b r \ 2n ti^_ H 

2n-2 



As b c ~ (s/M^ +2 ) 1//n , and b d = 1/Af s , we note that the cross section grows 
oc s 2 Mg n ~ 2 /M^, n+4 . Thus, although the cross section is comparatively small in this re- 
gion, it has a stronger growth rate oc s 2 than in region 1. 

For the differential cross section, we note that the t-channel Born amplitude is pro- 
portional to 1/A; 2 for —k 2 3> M 2 . This implies that the cross section has a forward peak. 
It corresponds to scattering at distances small compared to 1/M S , in the 1/r potential 
from eq. ( [4.6[ ). There is however no forward divergence since the growth is softened at 
—k 2 ~ M 2 , i.e. at distances comparable to the brane thickness. 



6.3 Region 3, s < M 2 and |A Bo m^| < 1 

In region 3 the cross section is also dominated by the Born amplitude. But in this case 
the scattering is almost isotropic (apart from spin dependences) as the factor —k 2 in the 
propagator is small compared to the heavier and most important KK modes. This implies 
that we may also have important contributions from u- and s-channel exchanges. For 
identical particles, the u-channel contribution has the same magnitude as that from t- 
channel. 
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6.4 Region 4, s < M 2 , |A Bor n^| > 1 and Im(X) > Re(X) 

The one- loop t-type contribution in fig. ||a, is dominated by the imaginary part, originating 
from real intermediate states. If loop diagrams of this type dominate, the all-loop amplitude 
is approximated by the geometric sum in eq. ( |5.4| ). As in region 3, the result is then 
approximately isotropic, but here higher order corrections give some suppression compared 
to the Born approximation. For identical particles the u-type ladder is identical to the 
t-type ladder and hence equally important. For particle-antiparticle scattering s-channel 
contributions have to be considered. 

6.5 Region 5, s < M 2 a , |A Bo mX| > 1 and Im(X) < Re(X) 

In this region the one-loop diagram has a dominant real part. This implies that virtual 
intermediate states and inelastic reactions are important. Therefore non-ladder diagrams 
are expected to give large contributions, and we showed two examples in fig. |6|. This 
region is consequently much more complicated than the other kinematical regions. It 
corresponds to situations where the effective cut-off M s is large ("narrow brane" or strong 
"brane tension") and the energy is in an intermediate range. From fig. |3] we see that for, 
e.g. n = 3, M s must be larger than 5 Mq. In this paper we will not make any specific 
predictions for what might be expected in this kinematic region. 

7. Conclusions 

In the ADD model it is assumed that standard model particles live on a 4-dimensional 
brane, embedded in a (4 + n)-dimensional space with n compactified dimensions. In these 
only the gravitational field is allowed to propagate. If the brane is infinitely thin and 
infinitely rigid, the exchange of very massive Kaluza-Klein modes represents a contact 
interaction of infinite strength between the standard model particles. This is not physically 
acceptable and different ideas have been proposed to regularize the scattering process. 

If the brane has a finite width, or if it is not infinitely well localized, the exchange of 
KK modes will be suppressed for KK wavelengths shorter than the width of the brane, or 
the size of its fluctuations. This will therefore give an effective cut-off (denoted M s ) for 
high KK masses, which does not have to be of the same magnitude as the fundamental 
Planck mass Mq. 

In this paper we have studied the effect of such a cut-off on the scattering of stan- 
dard model particles at various energies. We find that several troublesome infinities and 
divergencies are removed. The scattering process depends on three different energy scales, 
the collision energy y/s, the fundamental Planck scale Mo, and the cut-off scale M s . The 
Planck scale, Mb = (87rit m G , 7v(4) _1/ ^ n+2 ^, depends on the compactification radius R of 
the extra dimensions and the magnitude of Newton's constant, while the effective cut-off 
depends on the width of the brane, M s ~ (brane thickness) -1 , or the fluctuations in its 
position. These scales are thus not automatically related. Clearly the compactification 
scale R must be larger than the brane width 1/M S . 

Depending on the relative magnitude between these scales, we have here studied five 
different kinematical regions with different dynamical behavior. In one region (region 1 
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in fig. H), the scattering is dominated by small angles, and the eikonal approximation is 
applicable. Here we recognize classical scattering in a \/r n+l potential and the results of 
Giudice-Rattazzi- Wells [4]. In two other regions (2 and 3 in fig. ||) the Born approximation 
is applicable. In one of these (region 2) forward scattering dominates, and corresponds 
to scattering in a 1/r potential, but with a coupling enhanced by a factor proportional to 
^ compactmcation Radius ^ n CO mpared to scattering in the ordinary 1/r Newtonian large distance 
potential. In the other Born region (region 3) the scattering is approximately isotropic, 
as expected in [5,6]. In a fourth region the exponentiation from ladder-type diagrams 
in the eikonal region is replaced by a geometric sum. The scattering is expected to be 
mostly elastic since on-shell intermediate states dominate, but approximately isotropic. In 
the last region inelastic processes and non-ladder loop diagrams are important and make 
predictions very difficult. The boundaries between the different regions are expressed in 
the three mass scales involved, as illustrated in fig. ||. 
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